Limits on Topological Defects Models of the Ultrahigh Energy Cosmic Rays [] 

G\ ; 

^ \ R.J. Protheroe 1 and Todor Stanev 2 

■ 1 Department of Physics and Mathematical Physics, University of Adelaide 

' Adelaide, SA 5000, Australia 

, 2 Bartol Research Institute, University of Delaware, Newark, DE 19716, U.S.A. 

' Using the propagation of ultra high energy nucleons, photons and electrons in the universal radia- 

tion backgrounds, we obtain limits on the luminosity of topological defect scenarios for the origin of 
the highest energy cosmic rays. The limits are set as function of the mass of the X-particles emitted 
by the cosmic strings or other defects, the cosmological evolution of the topological defects, and the 

, strength of the extragalactic magnetic fields. The existing data on the cosmic ray spectrum and 

on the isotropic 100 MeV gamma-ray background limit significantly the parameter space in which 
topological defects can generate the flux of the highest energy cosmic rays, and rule out models with 
the standard X-particle mass of 10 16 GeV and higher. 
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^\ The cosmic ray events with the highest energies so far 
detected have energies of 2 x 10 11 GeV U and 3 x 10 11 

|H , GeV 0j . The question of the origin of these cosmic rays 

^ ■ having energy significantly above 10 11 GeV is compli- 
\ cated by propagation of such energetic particles through 

^ • the universe. The threshold for pion photoproduction on 
K^i | the microwave background is ~ 2 x 10 10 GeV, and at 
. 3 x 10 11 GeV the energy-loss distance is about 20 Mpc. 

C$ ' Propagation of cosmic rays over substantially larger dis- 
tances gives rise to a cut-off in the spectrum at ~ 10 11 
GeV as was first shown by Greisen || , and Zatsepin and 
Kuz'min f§, the "GZK cut-off". 

The standard cosmic ray acceleration mechanism, 
shock acceleration, leads to a power-law energy spec- 
trum, dn/dE oc E~ a , with differential index a > 2. To 
reach energies of ~ 10 11 GeV one requires the conditions 
present in powerful radio galaxies j^]. 

An alternative explanation of the highest energy cos- 
mic rays is the topological defect (TD) scenario 
where the observed cosmic rays are a result of top-down 
cascading, from the GUT scale energy of ~ 10 16 GeV or 
higher Go], down to 10 11 GeV and lower energies. Gen- 
erally, these models put out much of the energy in a very 
flat spectrum of photons and electrons extending up to 
the mass of the "X-particles" emitted. Approximating 
this spectrum by monoenergetic injection of photons of 
energy 10 15 GeV, Protheroe and Johnson JTTJ showed 
that spectra from single TD sources can not explain the 
2-3 x 10 11 GeV events. 

The main problem with topological defect models is 
the wide range of model parameters in which this sce- 
nario could, in principle, be applied. Parameters of TD 



scenarios include: mass of the X-particle, energy spectra 
and final state composition of the decay products, and 
cosmological evolution of the topological defect injection 
rate [ p2|Jl^ | . The problem of propagation is more severe 
than for the case of acceleration scenarios because most 
of the energy from X-particle decay emerges in electrons, 
photons and neutrinos, with only about 3% in nucleons. 
The electrons and photons initiate electromagnetic cas- 
cades in the extragalactic radiation fields and magnetic 
field, resulting in a complicated spectrum of electrons and 
photons which is very sensitive to the radiation and mag- 
netic environment. For example, recently the HEGRA 
group pi ] have placed an upper limit on the ratio of 7- 
rays to cosmic rays of ~ 10~ 2 at 10 5 GeV and, using a TD 
model calculation which neglected the IR background 
and gave a higher ratio, argued that TD models were 
ruled out. However, inclusion of the IR would reduce the 
10 5 GeV 7-ray intensity to well below the HEGRA limit. 

Protheroe and Johnson |l5| considered one set of pa- 
rameters (Mxc 2 — 10 15 GeV, constant injection per co- 
moving volume, B = 10~ 9 gauss) and ruled out TD as 
the origin of the 2-3 x 10 11 GeV events. This was mainly 
due to the high gamma-ray intensities at observable ener- 
gies in the electromagnetic cascade initiated by electrons 
and photons in the TD spectrum above 10 11 GeV. The 
unification mass obtained from an analysis of LEP data 
is io 16 0±0 - 3 GeV, and the X-particle mass cannot be 
far from this. For an X-particle mass close to the unifica- 
tion mass, i.e. higher than the 10 15 GeV used in ref. [flB] , 
even more energy would be injected into this cascade, 
and the gamma-ray intensities would violate the obser- 
vational contraints even more. Ref. |l5| has therefore 
already ruled out TD as the origin of the 2-3 x 10 11 GeV 
events. Recently, however, Lee [|l3| and Sigl, Lee and 
Coppi jl6| have claimed that lower X-particle masses are 
possible, and adopting Mxc 2 = 10 14 GeV, and a lower 
magnetic field, suggested the TD scenario is not ruled 
out. In this letter we consider several TD scenarios to put 
limits on the luminosity of the particle fluxes injected by 
topological defects as a function of the X-particle mass, 
the cosmological evolution of the topological defects and 
the strength of the extragalactic magnetic field, and con- 
sider for what range of parameters TD could explain the 
2-3 x 10 11 GeV events. We confirm the conclusion of 
Protheroe and Johnson (ll| that for X-particle masses 
of 10 15 GeV or higher TD can not explain the 2-3 x 10 11 
GeV events and severely limit models with lower Mx- 

We use the same injection spectra and TD evolution 
as in ref. fig ]. This is approximately an i? -15 spectrum 
extending up to ~ Mxc 2 /2 confining ~ 3% nucleons and 
97% pions. In the matter dominated era of the universe, 
and assuming go = 0.5, the injection rate per co-moving 
volume is Q(t) — Qo{t/t )~ 2+p where p = 1 for ordinary 



scribed in ref. Jll]]. The following processes are included: 
77 — > e + e _ on the microwave, radio and IR/optical back- 
ground, IC scattering on the same backgrounds, triplet 
pair production and double pair production on the mi- 
crowave background, synchrotron radiation in the extra- 
galactic magnetic field, and redshifting due to expansion 
of the universe. Nucleons undergo pion photoproduc- 
tion interactions and protons undergo Bethe-Heitler pair 
production in the same environment, and neutron pro- 
duction and decay are taken into account. 

For the radio background we use the spectrum of Clark 
et al. (IT]] . Other estimates of the radio background based 
on data on radio galaxies and ordinary galaxies give a ra- 
dio background extending to significantly lower frequen- 
cies , and we shall discuss the effect of using dif- 
ferent radio spectra elsewhere j2(j] . Magnetic field values 
we use are lCT 15 , 10" 12 , lCT 11 , . . . 1CT 8 G. The values at 
the high end of this range may be appropriate if topo- 
logical defects are seeds for the formation of galaxies and 
larger structures in the universe f2lf| where fields are gen- 
erally higher than average. For the infrared background 
we adopt a spectrum p2| which is based on the model 
of Stecker et al. [^3| but constrained at low frequencies 
by upper limits derived by us from the error bars on 
the microwave background measured by the FIRAS ex- 
periment on COBE At 3 x 10~ 3 eV, where the 
microwave background is decreasing rapidly with energy, 
our IR spectrum is a factor of 5 lower than that used by 
Lee fjj]. 

For a uniform distribution of topological defects we 
obtain the total intensity by integrating over redshift the 
results obtained for propagation over fixed distances, tak- 
ing account of topological defect evolution and cosmo- 
logical expansion assuming Ho = 75 km s _1 Mpc^ 1 and 
<7o = 0.5. The result for M x c 2 = 10 141 GeV, a magnetic 
field of 10~ 9 gauss, and p = 2, is shown in Fig. l] where 
we have normalized the spectrum of "observable parti- 
cles" (nucleons, photons, electrons) to the 3 x 10 11 GeV 
point (cosmic ray data are taken from p5| , and the high- 
est point is from Q). Lee jl3| has published a spectrum 
for similar input parameters and it is in acceptable agree- 
ment with the present work except for MeV-PeV 7-rays 
where our result is about a factor of 10 lower. We suspect 
this is because of our lower IR field, and this appears to 
be confirmed by results presented by Lee which show the 
7-ray intensity to be significantly lower if the IR field is 
neglected. With our lower IR field, and consequent lower 
7-ray intensity, we are less likely to rule out topological 
defect models due to excess 7-ray production. For nor- 
malization to the 3 x 10 11 GeV data, the injection rate of 
energy in X-particles would be ~ 6 x 10~ 35 erg cm -1 s -1 . 
Notice that above 10 11 GeV photons dominate the spec- 
tra of observable particles, and that over some ranges 



the topological defect models. Figure || shows the energy 
injection rate at the present epoch, such that the inten- 
sity of observable particles is normalized to the 3 x 10 11 
GeV point, as a function of Mx for various extragalactic 
magnetic fields and evolution models. 

Synchrotron radiation is very important in determin- 
ing the 7-ray spectrum at MeV-PeV energies which can 
vary by orders of magnitude depending on the magnetic 
field. Limits on the injection rate obtained from com- 
paring the predicted 0.1-10 GeV intensities with SAS-II 
|p6| and preliminary EGRET pT]] data are only lower 
than the injection rate obtained from normalization at 
3 x 10 11 GeV (and thus rule out a TD origin for the 
2-3 x 10 11 GeV events) for the highest magnetic fields. 
Where this limit is lower than the injection rate obtained 
from normalizing the intensity of observable particles to 
the 3 x 10 11 GeV point, these limits have been added to 
Figure § for the three evolution models. We see that the 
7-ray data provide the strongest constraint for models 
with high Mx, high B and weak evolution. No models 
with Mxc 2 < 10 144 GeV are excluded by the constraints 
imposed so far, so if we used only these two constraints 
we would agree with Sigl et al. |ll| that TD scenarios are 
not ruled out. 

A further constraint, not considered by Sigl et al. |l6| , 
comes from the intensity of potentially observable par- 
ticles above 3 x 10 11 GeV. This constraint has already 
been used by Protheroe and Johnson |D| to rule out the 
model with M x = 10 15 GeV, p = 2, and B = 10~ 9 gauss. 
Here we use the fact that 1 event was observed by the 
Fly's Eye between 10 1145 GeV and 10 1155 GeV, together 
with the published intensity at 10 115 GeV, to obtain the 
exposure factor of the Fly's Eye experiment at this en- 
ergy. Assuming the exposure factor has the same value 
also at higher energies (a reasonable assumption as at 
these energies optical transmission will limit the distance 
to observable air showers rather than the inverse-square 
law), we can estimate the number of events which should 
have been observed above 3 x 10 11 GeV. Given that no 
events have been seen above this energy, we set a 90% 
upper limit of 2.3 events which, when compared with the 
expected number of events, sets a new upper limit to the 
rate of injection of energy in X-particles. This limit is 
approximately independent of topological defect evolu- 



against Mx for various 
ris limit is lower than the 



tion and is plotted in Figure 
magnetic fields. In all cases t 
injection rate required to explain the 2-3 x 10 11 GeV 
events, and so it would appear that, subject to 7-rays 
and electrons above this energy being detectable by the 
Fly's Eye as discussed below, topological defect models 
are ruled out as the explanation of the 2-3 x 10 11 GeV 
events. Comparing Figs. 2 and 3, and extrapolating to 
10 16 GeV, it is obvious that TD models with standard 



trons which might be undetectable. Energetic 7-rays en- 
tering the atmosphere will be subject to the LPM effect 
p8[ (the suppression of electromagnetic cross-sections at 
high energy) which becomes very important. The ra- 
diation length changes as (E/E^pm) 1 ^ 2 , where E h p M = 
6.15 x 10 4 £ ra d GeV, and 4ad is the standard Bethe-Hcitler 
radiation length in cm 29 1. We find that average shower 
maximum will be reached below sea level for energies 
5 x 10 11 GeV, 8 x 10 11 GeV, and 1.3 x 10 12 GeV for 
gamma-rays entering the atmosphere at cos 9 = 1, 0.75, 
and 0.5 respectively. Such showers would be very diffi- 
cult to reconstruct by experiments such as Fly's Eye and 
at best would be assigned a lower energy. In this case, 
we should treat electrons and 7-rays as unobservable, and 
normalize the nucleon intensity to the 3 x 10 11 GeV data. 
This has the effect of increasing the predicted 7-ray in- 
tensities, and the new upper limits to the rate of injection 
of energy in X-particles would be as given in Figure ||. 
We now see that normalizing to the 3 x 10 11 GeV data 
violates the 7-ray data for all models with p = 2, models 
with p = 1 and M x > 10 131 -10 137 GeV depending on 
B, and models with p = and M x > 10 139 -10 14 - 9 GeV 
depending on B. Thus, models with standard Mx would 
also be ruled out as the explanation of the 2-3 x 10 11 GeV 
events. 

Before entering the Earth's atmosphere 7-rays and 
electrons are likely to interact on the geomagnetic field 
(see Erber [[30| for a review of the theoretical and ex- 
perimental understanding of the interactions). In such 
a case the 7-rays propagating perpendicular to the ge- 
omagnetic field lines would cascade in the geomagnetic 
field, i.e. pair production followed by synchrotron ra- 
diation. The cascade process would degrade the 7-ray 
energies to some extent (depending on pitch angle), and 
the atmospheric cascade would then be generated by a 
bunch of 7-rays of lower energy. Aharonian et al. ]3l| ] 
have considered this possibility and conclude that this 
bunch would appear as one air shower made up of the su- 
perposition of many sub-showers of lower energy where 
the LPM effect is negligible, the air shower having the 
energy of the initial 7-ray outside the geomagnetic field. 
If this is the case, then 7-rays above 3 x 10 11 GeV would 
be observable by Fly's Eye, etc., and the upper limits 
presented in Figure |3J would stand, ruling out a TD ori- 
gin for the 2-3 x 1IT 1 GeV events. There is however, 
some uncertainty as to whether pair production will take 
place in the geomagnetic field. This depends on whether 
the geomagnetic field spatial dimension is larger than the 
formation length of the electron pair, i.e. the length re- 
quired to achieve a separation between the two electrons 
which is greater than the classical radius of the electron. 
This question of whether or not pair production in the 
geomagnetic field takes place needs further investigation. 
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FIG. 1. Spectra at Earth for the topological defect model 
discussed in the text. SAS-2 and EGRET 7-ray data are 
shown at GeV energies, and HEGRA data at 100 TeV. 
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FIG. 2. Maximum rate of injection of energy in 

X-particles as a function of Mx for various magnetic fields 
and evolution models based on normalization of predicted in- 
tensity of "observable particles" to the 3 x 10 11 GeV point, 
or using the 7-ray data as upper limits (the lower of the two 
is plotted). Numbers attached to curves give log(B/ 1 G). 
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FIG. 3. Maximum rate of injection of energy in 

X-particles as a function of Mx for various magnetic fields 
and evolution models based on the non-observation of cosmic 
rays above 3 x 10 11 GeV. 
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FIG. 4. Maximum rate of injection of energy in 

X-particles as a function of Mx for various magnetic fields 
and evolution models based on normalization of predicted in- 
tensity of nucleons (solid curves), and using the 7-ray data 
as upper limits (dashed curves) for B = 1CT 15 G (highest 
curves) to B = 1CP 8 G (lowest curves). 
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